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Summary 

Often, the shear and tensile forces and moments are extracted from finite element analyses to be used 
in off-line calculations for evaluating the integrity of structural connections involving bolts, rivets, and 
welds. Usually the maximum forces and moments are desired for use in the calculations. In situations 
where there are numerous structural connections of interest for numerous load cases, the effort in finding 
the true maximum force and/or moment combinations among all fasteners and welds and load cases 
becomes difficult. The Element Load Data Processor (ELDAP) software described herein makes this 
effort manageable. This software eliminates the possibility of overlooking the worst-case forces and 
moments that could result in erroneous positive margins of safety and/or selecting inconsistent 
combinations of forces and moments resulting in false negative margins of safety. In addition to forces 
and moments, any scalar quantity output in a PATRAN (MSC Sofware Corporation) report file may be 
evaluated with this software. 

This software was originally written to fill an urgent need during the structural analysis of the Ares 
I-X Interstage segment. As such, this software was coded in a straightforward manner with no effort 
made to optimize or minimize code or to develop a graphical user interface. 

Symbols 

A percent yield at maximum preload for yield strength interaction; (percent of Ft y ) / Ft u for 

ultimate strength interaction. See Equations (19) to (21). 

At fastener tensile area 

B n<\)SF / Ft, At for yield strength interaction; n§SF / Ft u At for ultimate strength interaction. 

See Appendix A. 

C SF/0.571Ft y At for yield strength interaction; SF /0.577 Ft u At for ultimate strength 

interaction. See Appendix A 
F r resultant force 

F ty fastener material tensile yield stress 

Ft u fastener material ultimate tensile stress 

F x ,y,z force components 

[IM] matrix of stress ratios 

M r resultant momentum 

M x ,y,z moment components 

n loading plane factor (Ref. 1) 

R a .b,s stress ratios due to axial, bending, and shear loads, respectively 

’NASA Glenn Research Center’s Lewis’ Educational and Research Collaborative Internship Project (LERCIP). 
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Rx,y,z stress ratio components 

SF yield factor of safety for yield strength interaction; ultimate factor of safety for ultimate 

strength interaction (Ref. 1) 

[Shear] matrix of shear values 

V x ,y,z shear force components 

x, y, z force coordinates 

()) load factor (Ref. 1) 


Introduction 

Bolts are frequently modeled in NASTRAN (MSC Sofware Corporation) as rigid, beam, or spring 
finite elements for the purpose of extracting forces and moments to be used in subsequent off-line bolted- 
joint calculations; and also in the case of beam and spring elements, to model the bolt or bolted joint 
stiffness. Typically, the grid point (GP) force or multipoint constraint (MPC) force components are 
extracted at each GP associated with the rigid finite element representing the fastener. Element forces 
may be extracted from the beam and spring elements if they are used. Typically, the square root of the 
sum of the squares (i.e., the root-sum-square (RSS)) of two of the force components is computed to arrive 
at the shear force on the fastener. The third force component is usually the fastener axial force. Similar 
approaches may be used to evaluate welded and riveted connections. 

To obtain fastener forces, especially for multiple load cases, it is convenient to generate a report file 
using PATRAN (MSC Software Corporation), listing the GP or MPC forces as shown in Figure 1. For 
each load case, the report file can be made to list the x,y ^-coordinates of the GPs, as well as the GP x, y, 
and z force components. PATRAN automatically lists the GP number (denoted in the PATRAN report 
file as “Entity ID”). Concluding the report file is a summary of the maximum force components 
correlated with their corresponding GP and load case. The x, y, and z force components may be output 
with respect to any coordinate system within PATRAN. 

The Element Load Data Processor (ELDAP) was originally intended for use solely on GPs and 
consequently is limited to reading seven different data quantities, which were envisioned to cover GP 
identification number “ID,” three coordinates, and three forces or moments as previously mentioned. 
However, ELDAP can accommodate element forces as well. When element forces are requested in the 
PATRAN report file, PATRAN automatically outputs an element position ID along with the element ID, 
consequently taking up two data quantities and leaving five remaining data quantities choices. For 
example, when working with element forces the PATRAN report file could be specified to list any two of 
its coordinates (x, y, z), and its three x, y, and z force components. The element (Entity) ID and position 
ID are automatically generated, resulting in seven data quantities. 

Generally, the maximum force components do not result from the same load case or GP (fastener) 
location. This is shown in Figure 1, where load cases 17, 5, and 19 produce the maximum x, y, and z force 
components, respectively, at GPs 94467, 102339, and 136846, respectively. It is straightforward to use 
the maximum force components listed in the report summary to calculate the fastener shear force, even if 
the forces did not originate from the same load case and/or GP. Doing this results in what is known as a 
maximum-on-maximum load, often referred to as a “max -on-max load,” giving the largest shear force. 
This approach, being straightforward and convenient, will yield conservative bolted-joint margins of 
safety. However, this may result in the unnecessary selection of higher strength fasteners and possible 
false negative margins. 
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Figure 1. — PATRAN (MSC Software Corporation) report file. 
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Figure 1. — Continued. 
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Figure 1. — Concluded. 


A list of the load cases containing one of the maximum force components is provided in the 
PATRAN report summary. One should not usually RSS the two force components from a load case listed 
in the PATRAN report summary, as this generally will not yield the maximum shear force. As an 
example, consider the hypothetical load case summary as shown in Figure 2. Here, the load cases that 
contain one of the maximum force components are listed. There are four load cases total, meaning one 
load case is not listed because it does not have a maximum force component. Here the maximum forces 
(fastener loads) in the x-, y-, and z-directions occur in load cases 2, 3, and 1, respectively. Assume the 
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Nonlisted worst load case 

Load case 4: F x = 3 F y - 2.5 F z = 2.9 \/ x = 3.83 

Figure 2. — Hypothetical load summary where F xyz are force components 
and V x is shear force in y,z-plane. 

analyst wants to calculate the shear force in the y,z-plane. If the analyst performs an RSS of the y and z 
force components from the same load case, the resultant shear force from load cases 1, 2, and 3 would be 
3.16, 3.2, and 3.28 lbf, respectively. However, these load cases do not produce the largest shear force. 

This occurs in load case 4 where the shear force is 3.83 lbf. Load case 4 was not listed in the summary 
because its force components were not a maximum. Therefore, it is also possible to use loads that are too 
low resulting in false positive margins of safety. Note that the examples in this report use English units 
(lbf, in.) though metric units apply (N, m) and could be used as well. 

Element Load Data Processor (ELDAP) is a program based on MATLAB (The Mathworks, Inc.) that 
was developed to directly address these problems. The purpose of the program is to provide a tool that 
will enable the analyst to determine the maximum shear load and corresponding axial load from among 
multiple load cases and grid points using the PATRAN-generated report file. Specific GPs (fasteners) 
may be selected, shear forces in different planes and their corresponding axial forces may be obtained, 
and GP (fastener) loads and positions may be plotted. This program will enable the analyst to compute 
more-accurate bolted-joint margins by avoiding possible false negative margins as a result of using max- 
on-max loads, and by avoiding missing the true maximum bolt loads arising out of the use of the 
summary load cases containing one maximum force component. 

Equations 

Originally the program was developed to calculate fastener shear forces (by use of an RSS of two scalar 
forces) and also provide the third consistent scalar quantity, fastener axial load, associated with the two RSS 
components; consequently this report will focus on these aspects as the primary example. However, the 
program can calculate the RSS of any two or three scalar quantities that PATRAN provides as output. 

This section will present the equations used in this program to calculate the forces and moments from 
the ELDAP input file (PATRAN report file). As previously mentioned, fastener shear forces will be the 
focus of the example. However, as shown below, the equations are applicable to moments as well as 
forces. ELDAP first reads the PATRAN report file and converts its data into matrix form to be used in 
subsequent calculations as shown in Equations (3) to (13). Equations (3) to (5) are used to calculate the 
shear forces V x ,y,z from the three force components F x , F y , and F z . Subscripts on V refer to the axis normal 
to the plane of shear (refer to Fig. 3). Equation (6) calculates the resultant of all three force components. 
By replacing F x , F y , and F z with M x , My, and M z , respectively, the moment components may be RSS’d as 
shown in Equations (7) to (10). 
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Up to this point, fastener axial and shear loads have been treated separately and not in combination. 
The interaction equation in Reference 1 combines stress ratios R due to axial, bending, and shear loads: 

(R a +Rbf + R]< 1 (1) 

ELDAP can determine the worst-case fastener loads based on a combination of both axial and shear 
loads, also represented by an interaction equation. The interaction Equations (1 1) to (13) are derived in 
Appendix A and are based on the equations in Reference 1 . The interaction equations of Reference 1 
include the effects of bolt bending; however, the interaction Equations (1 1) to (13) used in the present 
work omit the bolt bending term for simplicity. 


R} + R} < 1 ( 2 ) 

In accordance with Reference 1 it is desired that the quantities R x , R y , and R z be less than or equal to 1 
for fastener shanks parallel to the x-, y-, or z-directions, respectively. 

The constants A, B, and C are functions of the material properties and geometry of both the fastener 
and the joint and are defined in the Nomenclature section of this report and derived in Appendix A. These 
three constants are input into ELDAP if the calculation of the interaction equations is desired. 


V x =p 2 +F 2 

(3) 

V y =^F 2 + F 2 

(4) 

Vz=^F} + F 2 

(5) 

F r =p 2 + F 2 +F 2 

(6) 

M yz =^M 2 +M 2 

(7) 

m xz =^m 2 +m} 

(8) 

+ 

ll 

5 

(9) 

Mr = ^M 2 +M 2 +M} 

(10) 

R x ={A + BF X ) 2 + \fp 2 + F 2 J 

(11) 

Ry={A + BFy) 2 + \c4f 2 + F 2 )' 

(12) 

R z = (A + BF Z ) 2 + \Cyj F} + F 2 f 

(13) 
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The quantities calculated in Equations (3) to (6) and (7) to (10) are shown in Figures 3 and 4, 
respectively. Equations (3) to (6) are used to form four individual vectors in ELDAP: V x , V y , V z , and F,. 
which are subsequently arranged in columns to form the matrix of shear values: [Shear] = [V x V y V z F r \, as 
shown in Figure 5. If moment quantities were generated in the PATRAN report file, this matrix would 
consist of four column vectors of moments \M yz M xz M X y M r \. 

ELDAP determines the maximum value of Equations (3) to (13) and reports them in 1 1 output files, 1 
each for the 1 1 calculated quantities. The files are described in the section “Output Files.” 




/ 


Figure 3. — Shear forces V xyz and resultant 
force F r with respect to force components 
F x,y,z- 


z, M z 


x, 





Figure 4. — Moments M xyxzyz and M r with 
respect to moment components M x y z . 
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These are the calculations made from the last three 
columns of each load case (matrix M) 


Vx 


1336.910700 

1336.910532 

921.978301 

921.977126 

2996.716673 

2996.716677 

1569.413657 

1569.414514 

4438.840459 

4438.840479 

1565.275950 

852.244167 

1605.490155 

1138.256118 

1422.372006 

1472.718112 

5145.043366 

668.093756 

668.097576 

1565.277390 

852.247083 

1605.491887 

1138.258460 

1422.374163 

1472.720095 

5145.043342 

5566.308307 

5566.308355 

5573.879569 

5298.053740 

520.352369 

5573.879614 

5745.425428 

5298.053742 

520.356588 

5745.425433 

5288.944781 

3841.968194 

5904.629305 

5288.944799 

3841.968245 

5904.629315 

5198.865162 

457.163287 

5766.569289 

5198.865148 

2174.805757 

457.167542 

5766.569288 


Vy 

1279.983308 
1279.983483 
449.532650 
449.535043 
2993.515464 
2993.515460 
1274.295747 
1274.294678 
4379.917682 
4379.917661 
1357.112018 
769.407175 
1426.171544 
1057.540873 
1304.158113 
1321.714669 
5040.770290 
548 .199774 
548.195122 
1357.110349 
769.403948 
1426.169579 
1057.538351 
1304.155762 
1321.712455 
5040.770303 
5457.745959 
5457.745930 
5479.005194 
5296.683299 
379.036501 
5479.005153 
5739.797243 
5296.683297 
379.030721 
5739.797239 
5230.297284 
3839.214855 
5903.688393 
5230.297268 
3839.214804 
5903.688382 
5166.131962 
324.429306 

5766.506567 
5166.131974 
2157.782563 
324.423308 

5766.506568 


Vz 

386.949527 
386.949524 
807.199307 
807.199298 
138.565387 
138.565387 
918.530960 
918.530940 
720.894669 
720.894666 
790.693334 
455.491546 
756.400025 
495.348334 
613.808841 
679.184347 
1030.597532 
467.235436 
467.235441 
790.693320 
455.491551 
756.399998 
495.348330 
613.808843 
679.184339 
1030.597470 
1094.117333 
1094.117430 
1024.956806 
120.642423 
434.412771 
1024.956830 
254.394740 
120.642421 
434 . 412781 
254.394743 
787.965145 
161 . 176558 
123.762789 
787.965162 
161.176554 
123.762793 
586.105688 
401.946911 
43.066152 
586.105675 
271 . 924309 
401.946909 
43.066150 


Fr 

1337.052248 
1337.052247 
922.955663 
922.955655 
2996.718707 
2996.718707 
1570.125473 
1570.125462 
4438.844030 
4438.844030 
1567 . 966482 
873.435764 
1609.923120 
1153.123944 

1431.909879 
1479.378194 
5145.044589 

694 . 687416 
694.687419 
1567.966475 
873.435767 
1609.923108 
1153.123942 

1431.909880 
1479.378190 
5145.044576 
5566.321580 
5566.321599 
5573.964829 
5298.055398 

549.158316 

5573.964833 

5745.428723 

5298.055398 

549.158324 

5745.428723 

5289.132048 

3842.282407 

5904.807410 

5289.132050 

3842.282407 

5904.807410 

5199.069065 

487.757094 

5766.618335 

5199.069064 

2174.827407 

487.757092 

5766.618335 


Figure 5. — ELDAP output file of matrix of shear values, [Shear], 


User’s Section 

ELDAP is an interactive program running in MATLAB. It requires a PATRAN report file and 
optionally a text file of GPs arranged in one row and delimited by commas. This section will present the 
steps necessary to generate the ELDAP input files. The program input prompts and graphing options will 
be discussed. An example of the required input file as well as the output files will be shown. In this report 
“node” will be synonymous with “grid point.” 
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Generating the PATRAN Results (ELDAP Input) File 

The analyst must request the appropriate output quantities (GP forces, MPC forces, or element forces) 
in the NASTRAN bulk data file to be able to capture them in a PATRAN report file. This section will 
demonstrate the on-screen process for calculating fastener forces using GP forces and MPC forces. Note 
that terms in brackets ( ) are buttons, bold terms are categories, and italic terms are options to select. 

Place the desired GPs and finite elements representing the fasteners in the PATRAN display window. 
This may be accomplished by creating a group within PATRAN consisting solely of the finite elements 
and their GPs representing the fasteners. If GP or MPC forces were requested as output quantities in 
NASTRAN, all that is needed in the PATRAN group are the GPs. Shown in Figure 6 is an example of 
GPs representing a pattern of bolts. In PATRAN, select the (Results) button and choose Action: Create, 
Object: Report, and Method: Overwrite File as shown in Figure 7. Continuing to refer to Figure 7, select 
the desired result cases in the Select Result Cases(s) list box. Select the desired output quantities: Grid 
Point Forces, MPC Forces from the Select Report Result list box. Select the quantities: X Location, Y 
Location, Z Location, X Component, Y Component, and Z Component from the Select Quantities list box. 



Figure 6. — PATRAN (MSC Software Corporation) viewport showing rigid (RBE2) elements used to model fasteners. 
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Action: Create 


Object: Report ^ | 

Method: Overwrite File w \ 

m ^ ^ 


Select Result Case(s) 



Select Report Resui 
Displacements. Translational 
Grid Point Forces, Applied Loads 
Grid Point Forces, Constraint Forces 
Grid Point Forces, Internal Forces 


S 3 



Appty | Reset 


| Group | ~ 


Select load cases 


Select: 

Grid point forces 
MPC forces 


Select: 

X Location 

Y Location 
Z Location 

X Component 

Y Component 
Z Component 


Figure 7. — Processing PATRAN (MSC Software 
Corporation) results to generate ELDAP input file. 
Selecting case(s), report result, and quantities. 


As shown in Figure 8, select the (Target Entities) icon and select Current Viewport from the drop- 
down box. Under Addtl. Display Control, select Nodes from the drop-down box. 

Referring to Figure 9, select the (Display Attributes) button and choose (File. . .). If necessary, create 
the folder where the report file will reside and type in the desired name of the report file. For Report 
Type, select Full from the drop-down box, which will generate all of the GP/MPC forces for all of the 
GPs residing in the viewport window or current group for the selected load cases. Select (Apply) to 
generate the report file. Note that the forces will be generated for all GPs in the active group even if the 
analyst has zoomed in and the GPs are not all showing in the viewport window. 

Figure 1 is an excerpt from a PATRAN-generated results file, which is the EFDAP input file. Most of 
the file has been eliminated for clarity. The first several lines of data for selected load cases are shown, 
followed by the Summary Information section containing the maximum force components and their 
corresponding load case. Fastly, an index correlating load case number with load case name is presented, 
most of which has also been eliminated for clarity. 
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Action; Create w | 

Object: Report w | 

Method: Overwrite File 


mM & 


(Target entities) icon — 

Select current viewport 
and nodes to capture 
rigid (RBE2) forces — cl' 


Target Entity 
Current Viewport 


Addtl. Display Controt 
Nodes w | 


Resuts I Group 


|p«jp I 


Figure 8. — Processing PATRAN (MSC Software Corporation) 
results to generate ELDAP input file. Selecting target entities. 


(Display attributes) icon 

Select Full to generate 
file containing all forces 

Apply | Reset 


Resets | Group | 



Action: Create w \ 

Object: Report » | 

Method: OverwrteFte ^ 


MjiM SB 


| C 'eackuptlnterstage.AdaptertlSI _C0f 

Format... 

Sorting Options 


Report Type: FuH w | 


Figure 9. — Processing PATRAN (MSC Software Corporation) 
results to generate ELDAP input file. Selecting report type. 
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Prompts 


This section is a brief overview of the input prompts that appear on the MATLAB window when 
ELDAP is run. As previously mentioned, this discussion will focus on extracting and calculating fastener 
loads, but the program may be used for any other available PATRAN output quantity (whether or not it 
makes sense). Figure 10 shows the MATLAB command line prompts and entries for this program. The 
first of the prompts will ask the user to enter the name of the file that contains the PATRAN report 
containing the GP, MPC, or element force data, followed by prompts for the number of load cases, and 
the number of nodes (GPs). The next prompt will ask the user if it is desired to calculate the fastener shear 
and tension interaction equation for each fastener. If so, the user is next prompted for the parameters A, 5, 
and C as defined in the Symbols list at the beginning of this report. The user can enter the values of A, B, 
and C based on either yield or ultimate strengths (even though the prompt for A refers to yield). ELDAP 
is capable of processing any number of load cases and nodes (fasteners), subject to the limitations of 
NASTRAN, PATRAN, and MATLAB. As presented in the section “Output Files,” ELDAP determines 
the maximum values of V x , V y , V z , F r (or M yz , M xz , M xy , Mr), R x , R y , and R z from among all load cases and 
GPs using Equations (3) to (13) and reports them along with their other corresponding consistent forces 
or moments in output files, one each for the calculated quantities. Additionally, ELDAP outputs the 
matrix of shear values, [Shear] = [V x V y V z F r \, and the matrix of stress ratios, [ IM\ — \R x R y R z ]. If 
moments are being evaluated, the matrix of shear values contains the moments, [Shear] = \M xy M xz M yz 
M r \. The subsequent prompts give the user the option of adding a standard prefix to the default file names 
(the default file names would then be a suffix) and asks the user to enter the desired prefix or hit return to 
accept the default file name listed on the screen. The default file names are Vcalc.out, Vx.out, Vy.out, 
Vz.out, Fr.out, Fx.out, Fy.out, Fz.out, Rcalc.out, Rx.out, Ry.out, and Rz.out. The user also has the option 
of changing the default name of each separate output file to a chosen name (e.g., Myname.out). If the 
prefix option is chosen, it will be used either for the default names (e.g., prefixVx. out) or the chosen 
names (e.g., prefix Myname.out). 

The subsequent prompts ask the user to enter the name of a file containing a list of specific GPs to be 
evaluated and if there is not one, to hit return. If a file name is entered, ELDAP makes calculations only 
on the GPs contained in that file. To use this option, the analyst must first construct a text file by creating 
a list of the desired GPs, making sure to delimit the list with commas and limit the list to one row of text. 

Next, ELDAP prints to the screen the maximum values and corresponding load cases for V x , V y , V z , 

F r , R x , R v , R z , F x, F y , and F z . 

The remaining prompts allow the analyst to make three types of plots of the processed data. The first 
of these prompts enables the analyst to make a three-dimensional (3D) plot of the original and processed 
data. To do this the analyst is asked to enter the load case from which to plot the data and also to enter the 
specific quantities to plot on the x-, y-, and z-axes. The user can refer to the data just printed to the screen 
to help choose which load case is of most interest. When entering the quantities to be plotted, the analyst 
needs to be aware of the orientation of the model to be sure to plot what is desired. ELDAP gives the 
analyst the option of making up to 15 consecutive plots. The second plot request will produce a graph of 
the location of the GPs along with their corresponding GP number. The third plot request will provide a 
graph of any one of the calculated quantities from Equations (3) to (13) for a particular GP (fastener) for 
each load case. 

The last ELDAP request will provide the opportunity to print all of the output files to the screen in 
concatenated form. 
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Please input the name of the file from Patran to be opened (extension required) ELDAPInput . rpt 
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Figure 10. — MATLAB (The Mathworks, Inc.) screen capture of ELDAP input prompts. Note that entries typed in by 
user are in italic font, and (Enter) indicates user pressed “Enter” button (this is “Return” on some computers) 
to select default entry. 
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Figure 10. — Continued. 
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Figure 10. — Continued. 
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Figure 10. — Concluded. 

Output files 

ELDAP determines the maximum values of V x , V y , V z , F r (or M yz , M xz , Mxy, M r ), R x , Ry, and R z from 
among all load cases and GPs using Equations (3) to (13) and reports them along with their other 
corresponding consistent forces (or moments) in output files one each for the calculated quantities. The 
corresponding default file names are Vx.out, Vy.out, Vz.out, Fr.out, Fx.out, Fy.out, Fz.out, Rx.out, Ry.out, 
and Rz.out. If moments instead of forces are being evaluated, the output files Vx.out, Vy.out, Vz.out, and 
Fr.out contain moments instead of forces. The first output file that EFDAP generates has the default name of 
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Vcalc.out and contains the matrix of shear values: [Shear] = [ V x V y V z F r \. This four-column matrix consists 
of the shear forces V x , V y , and V z in the y,z-, x,z-, and x,y-planes, respectively, as well as F r , which is an 
RSS of the three force components, as shown in Equations (3) to (6). If moments are being evaluated, this 
matrix contains the moments as calculated in Equations (7) to (10). Figure 5 is an example output file 
containing the matrix of shear values. The file has been shortened to fit on the page. ELDAP also calculates 
the interaction equation stress ratios R x , R y , and R z as shown in Equations (1 1) to (13). These stress ratios are 
stored in the three-column matrix [IM] = [R x R y R z ], This matrix is also printed to an output file with a 
default file name of Rcalc.out and in the same format as the matrix of shear values. The user has the option 
of naming these output files (e.g., MyName.out) or can accept the default name (e.g., Vx.out) given by the 
program. As previously mentioned, the user may also specify a prefix to the chosen or default file name. 

The default file names are included in the prompts as shown in Figure 10. 

Figure 1 1 is an example output file for V x , the shear force in the y,z-plane. Output files for V y , V z , F r , 
F x , F y , F z , Ry, and R z have identical formats. In the upper left of Figure 1 1 is the PATRAN report file load 
case number and row number for the data in this load case corresponding to the maximum shear force V x . 
Below that is the maximum shear force value. The NASTRAN load case that produced the maximum V x 
is provided next in the center of the figure. After that is a reprint of the specific row of data from the 
PATRAN report file (containing the maximum V x shear force) along with the RSS values of the 
combinations of the x-force (XComponent), y-force (YComponent), and z-force (ZComponent) 
components, denoted as V x , V y , V z , and F r per Equations (3) to (6). If moment results were used instead, 
the output file would have the moments about the x-, y-, and z-axes replacing the force components, and 
moment RSS values per Equations (7) to (10) replacing the RSS of the forces. Figure 1 1 presents the 
entire output file. 

Similarly, an output file for F x , the maximum axial force on the fastener, is also to be generated. The 
file also includes the same type of data and has the same format as that shown in Figure 1 1 . The specific 
row of data from the PATRAN report file containing the maximum force F x , along with the RSS values 
of the three combinations of the x-force (XComponent), y-force (YComponent), and z-force 
(ZComponent) components denoted as V x , V y , and V z is also provided, where the shear force V x 
corresponds to and is consistent with the axial force F x . In this way both types of output files, one for V x 
and one for F x , provide fastener maximum shear and correlated axial load and maximum axial load and 
correlated shear loads, respectively. 


Graphs 

ELDAP provides the analyst several plotting options. The first of these options is to make a 3D graph 
of any combination of the 10 data quantities (x, y, z, V x , V y , V z , F r , F x , F y , and F z ) from a particular load 
case previously generated. Figure 12 is an example 3D plot showing fastener (GP) shear forces (vertical 
axis) at each fastener location (horizontal axes). 

Figure 13 is an example of the second plotting option. This option plots the location of each GP and 
displays the GP number beside its corresponding location on the graph. When the GP numbers displayed 
on the plot overlap, the analyst can select (click) and drag the numbers apart. 

The last plot option is shown in Figure 14. Here the analyst can choose to plot any of the 10 data 
quantities (x, y, z, V x , V y , V z , F r , F x , F y , and F z ) for any one node (GP) as a function of all of the load 
cases. In this example, the GP is listed at the top of the graph, the load is indicated on the ordinate, and 
load case number on the abscissa. 
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Figure 1 1 . — Example of ELDAP output file Vx.out, for shear force in y,z-plane, V x . 
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z-location, in. Shear force, 


Load case number 15 



Figure 12. — ELDAP three-dimensional plot showing fastener shear force V x at each fastener location. 
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Figure 14. — ELDAP shear load (Ibf) for node (grid point) 132049 as function of all load cases. 
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Conclusion 

A MATLAB (The Mathworks, Inc.) computer program has been developed to facilitate the 
determination of the maximum NASTRAN (MSC Software Corporation) grid point (GP) forces and 
moments from among numerous GPs and numerous loads cases. This software denoted as ELDAP 
(Element Load Data Processor) determines the correlated shear and tension forces for each GP. The 
maximum shear force and corresponding tensile force or the maximum tensile force and the 
corresponding shear force for each GP are determined. This eliminates the possibility of recovering 
overly conservative loads as would be the case when only retrieving and combining the maximum force 
and moment components irrespective of the associated GP and load cases. It also eliminates the 
possibility of false positive margins that result from overlooking the root-sum-square (RSS) of less than 
maximum force components where the RSS would result in maximum forces. As such, this software 
eliminates erroneous margins of safety resulting from the use of inappropriate forces and moments. This 
software may also be used to evaluate any scalar quantity output from NASTRAN and PATRAN (MSC 
Software Corporation), in addition to the GP forces and moments. 
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Appendix A. — Derivation of Interaction 
Equations Based on Fastener Properties 

When a fastener experiences multiple loads simultaneously, such as axial, bending, and shear loads, 
the combined effect of all of the loads on the fastener integrity can be determined using an interaction 
equation. A typical interaction equation for axial and shear loads, omitting bending for simplicity, is 
considered in this work and is shown in Equation (2): 

rZ+rH i (2) 


Here, R a is the ratio of the bolt axial load to the allowable axial load, and R s is the ratio of the bolt shear 
load to the allowable bolt shear load. The stress ratios can be written in terms of either yield or ultimate 
strength. The following derivation will be based on yield strength. Similar equations may be derived 
based on ultimate strength by replacing Ft y with Ft,, (except that PLDmax remains in terms of Ft, ) and 
by replacing the yield factor of safety with the ultimate factor of safety. The allowable axial yield load is 
the product of the bolt allowable yield stress Ft, , and the bolt tensile area At. 

The ratio of the bolt load to the bolt allowable yield load may be expressed as 


R 


a 


Pb 

Ffy At 


(14) 


where Pb is the bolt axial load due to both the maximum preload PLDmax, and external load P, 

Pb — PLDmax + n§(SF xp) 


(15) 


and where n, (j), and SF are the loading-plane factor, load factor, and yield factor of safety, respectively, 
which are defined in Reference 1. Substituting Equation (15) into (14), 


_ PLDmax n<\>(SF xp) 
Fty At Fty At 


Assuming the bolt shank is aligned with the z-axis, the external axial bolt load may be expressed as 


P = F Z 

The maximum bolt preload is often given as a percentage of allowable material yield stress: 

PLDmax = (percent of F ty )A t 

Substituting Equations (17) and (18) into Equation (16), the axial stress ratio becomes 

percent of F ty A t n<\>(SF x F z ) 

Ra — 1 — 

Fty At Fty At 

Recognizing the constant terms, Equation (19) can be simplified as 

R a = A + Bx F z 


(17) 


(18) 


(19) 


(20) 
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where the constant coefficients are 


A = percent of fastener yield strength at maximum preload 

Fly At 


The ratio of external shear force to the allowable shear force may be expressed as 


Rs 


SFxV 


0.577 F ty At 


(21a) 

(21b) 


( 22 ) 


where Vis the external fastener shear force and 0.577 Ft y At is the allowable shear force, where the 
allowable shear stress is given as 0.577 F ty . The factor 0.577 is an average of the allowable shear stress 
for alloy and carbon steels (0.6 Ft y ) and stainless steels (0.55 Ft x ) as given in Reference 2. Since the 

fastener shank is aligned with the z-axis as previously mentioned, the fastener shear force may be 
computed as the root-sum-square of the two force components in the x,y-plane as shown in Equation (23): 


V = yj Fx + Fy 

Substituting Equation (23) into (22), 

SF yj F} + F} 
0.577 F ty At 


(23) 


(24) 


Recognizing the constant terms and simplifying, Equation (24) becomes 

R s = C^F} + Fj (25) 

where 


C = 


SF 

0.577 Ft y At 


(26) 


Substituting Equations (20) and (25) into Equation (2), the axial and shear interaction equation for a 
fastener with its shank parallel to the z-axis may be expressed as 

R z = (A + B x F z f + (c^+F^J (27) 

Similar equations may be developed for fasteners with shanks parallel to the other two orthogonal 
directions as shown in Equations (11) and (12) of this report. 
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Appendix B. — ELDAP Program Source Code 

This appendix provides the source code for the Element Load Data Processor (ELDAP). This source 
code is written in the MATLAB programming language and was originally written to fill an urgent need 
during the structural analysis of the Ares I-X interstage segment. As such, this software is coded in a 
straightforward manner with no effort made to optimize or minimize code or to develop a graphical user 
interface. 
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(969) filel3 = assigned to the function opening the output file to print the max calcu 

for (Ry) 

(991) filel4 = assigned to the function opening the output file to print the max calcu 

for (Rz) 

(902) filel5 = assigned to the function opening the output file to print the IM Matrix 
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points to be calculated 

(1280) Hby = converting Kbi2 from a number to a string 

(306) L = number of load cases 

(857) Label = assigned to the row of labels from the original file 
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This matrix contains the results from the interaction calculations from the 


last three columns of each load case (which are contained in matrix 
The three columns of -IM- are Rx, Ry, and Rz 
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disp('The program did not recognize the input, 
disp ('Enter y or n in a lowercase letter.'); 
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VZPRE2= [PRE2 , ' Vz .out ' ] ; 

disp(' Enter a name for the file containing calculations for Vz. 
disp ( [ ' Default VZPRE2 
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FZPRE2= [PRE2 , 'Fz.out' ] ; 

disp( 'Enter a name for the file containing calculations for Fz 
disp ( [ 'Default ' , ' ( ' ,FZPRE2, ')']); 
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disp(' Enter a name for the file containing calculations for Vz 
VZMCAL=input (' Default ( Vz . out ) ' , ' s ' ) ; 

if isempty (VZMCAL) ; 

VZMCAL = 'Vz .out ' ; 
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disp ('Enter a name for the file containing calculations for Fx. 
FXMCAL=input (' Default (Fx.out) 
if isempty (FXMCAL) ; 

FXMCAL = 'Fx.out ' ; 
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disp(' Please wait while ELDAP makes Calculations 

disp ( ' ' ) ; 
disp ( ' ' ) ; 
disp 

disp ( ' ' ) ; 
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fprintf ( f ile3 , ' \nRow Data Containing the max value of Vx ' ) ; 
fprintf (file3, ' \n%s ' , Label) ; 

fprintf (file3, ' \n%8f %16f %14f %14f %14f %14f %14f %14f %14f %14f 
, [Matrix (MaxVx, : , LoadMaxVx) Shear (MaxVx, : , LoadMaxVx) ] ) ; 
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for p2=l : (FN2 ) ; 

TempRy=max ( IM ( 1 : p2 , 2 , n2 ) ) ; 
if TempRy>TemplO ; 
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fprintf ( f ilel2 , ' \nRow Data Containing the max value of Rx ' ) ; 
fprintf (filel2, ' \n%s ' , Labe2 ) ; 

fprintf (filel2, ' \n%8f %16f %14f %14f %14f %14f %14f %14f %14f %14f 
, [Matrix (MaxRx, : , LoadMaxRx) IM (MaxRx, : , LoadMaxRx) ] ) ; 
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fprintf (filel4, ' \nThe max value for Rz ' ) ; 
fprintf(filel4, ' \n%f ' , IM (MaxRz , 3 , LoadMaxRz ) ) ; 


fprintf ( ' \nThe max value for Rz ' ) ; 
fprintf ( ' \n%f ' , IM (MaxRz, 3, LoadMaxRz) ) ; 
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MaxFy=f lat2 ; 
LoadMaxFy=covel ; 
Temp7=TempFy ; 
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fprintf ( file7 \nRow Data Containing the max value of Fx ' ) ; 
fprintf (file7, ' \n%s ' , Label) ; 

fprintf (file7, ' \n%8f %16f %14f %14f %14f %14f %14f %14f %14f %14f %14f 
, [Matrix (MaxFx, : , LoadMaxFx) Shear (MaxFx, : , LoadMaxFx) ] ) ; 


f ile8=f open (FYMCAL, ' wt ' ) ; 

fprintf ( f ile8 , 'Load Case for the max of Fy 
fprintf (file8, ' \n%f ' , LoadMaxFy ) ; 
fprintf (' \nLoad Case for the max of Fy ' ) ; 
fprintf ( ' \n%f ' , LoadMaxFy) ; 
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fprintf ( file9 \nThe max value for Fz ' ) ; 
fprintf (file9, ' \n%f ' , Matrix (MaxFz , 7 , LoadMaxFz ) ) ; 
fprintf (' \nThe max value for Fz'); 
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disp ( ' ' ) ; 

disp('The program did not recognize the input, 
disp ('Enter y or n in a lowercase letter.'); 
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disp ( ' ' ) ; 

file4=fopen ( VYMCAL, ' wt ' ) ; 

fprintf (1, 'Load Case for the max of Vy ' 

fprintf ( 1 , ' \n%f ' , LoadMaxVy ) ; 
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fprintf ( 1 \nRow Data Containing the max value of Vz ' ) ; 
fprintf ( 1 , ' \n%s ' , Label ) ; 

fprintf (1, ' \n%8f %16f %14f %14f %14f %14f %14f %14f %14f %14f %14f 
, [Matrix (MaxVz , : , LoadMaxVz ) Shear (MaxVz , : , LoadMaxVz ) ] ) ; 
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fprintf ( 1 , ' \n%s ' , Label ) ; 

fprintf (1, ' \n%8f %16f %14f %14f %14f %14f %14f %14f %14f %14f %14f 
, [Matrix (MaxFz , : , LoadMaxFz ) Shear (MaxFz , : , LoadMaxFz ) ] ) ; 
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